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Major advances have been made in the field of gamma-ray bursts in the last two years. The successful discovery 
of X-ray, optical and radio afterglows, which were predicted by theory, has made possible the identification 
of host galaxies at cosmological distances. The energy release inferred in these outbursts place them among 
the most energetic and violent events in the Universe. Current models envisage this to be the outcome of a 
cataclysmic event leading to a relativistically expanding fireball, in which particles are accelerated at shocks and 
produce nonthermal radiation. The substantial agreement between observations and the theoretical predictions 
of the standard fireball shock model provide confirmation of the basic aspects of this scenario. The continued 
observations show a diversity of behavior, providing valuable constraints for more detailed, post-standard models 
which incorporate more realistic physical features. Crucial questions being now addressed are the beaming at 
different energies and its implications for the energetics, the time structure of the afterglow, its dependence on the 
central engine or progenitor system behavior, and the role of the environment on the evolution of the afterglow. 



1. INTRODUCTION 

The discovery of gamma-ray bursts (GRB) by 
the Vela military test-ban treaty satellites was 
announced in 1973 |45| , and was quickly con- 
firmed by Soviet Konus satellite measurements 
|52| . Then, for 23 years GRB remained essen- 
tially just that: brief outbursts of gamma-rays 
which pierced, for a brief instant, an otherwise 
pitch-black gamma-ray sky. An intense debate 
festered for a long time on whether they were ob- 
jects in our galaxy or at cosmological distances. 
The first major breakthrough came in 1992 with 
the launch of the Compton Gamma-Ray Obser- 
vatory, whose superb results are summarized in 
a review by Fishman & Meegan ||2^ . In particu- 
lar the all-sky survey from the BATSE instrument 
showed that bursts were isotropically distributed, 
strongly suggesting either a cosmological or an 
extended galactic halo distribution, with essen- 
tially zero dipole and quadrupole components. 
The spectra are definitely non-thermal, typically 
fitted in the MeV range by broken power-laws 
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whose energy per decade vFi^ peak is in the range 
50-500 KeV the power law sometimes ex- 
tending to GeV energies |Q. GRB appeared to 
leave no detectable traces at other wavelengths. 



except in some cases briefly in X-rays | 119| , 13[ . 
The gamma-ray durations range from lO^"* s to 
about 10'^ s, with a roughly bimodal distribu- 
tion of long bursts of tb>,2 s and short bursts of 
tb<28 |4|], and substructure sometimes down to 
milliseconds. The gamma-ray light curves range 
from smooth, fast-rise and quasi-exponential de- 
cay (FREDs), through curves with several peaks, 
to highly variable curves with many peaks |^,^ . 
The pulse distribution is complex 72|, and the 
time histories can provide clues for the geometry 
of the emitting regions [ pTl|2^ . 

Theoretically, it was clear from early on that 
if GRB are cosmological, enormous energies are 
liberated in a small volume in a very short time, 
and an e^— 7 fireball must form ^8 32 ,117], which 
would expand relativistically. The main difficulty 
with this was that a smoothly expanding fire- 
ball would convert most of its energy into kinetic 
energy of accelerated baryons (rather than pho- 
tons), and would produce a quasi-thermal spec- 
trum, while the typical timescales would not ex- 
plain events much longer than milliseconds. This 
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problem was solved with the introduction of the 
"fireball shock model" |98|j57| ], based on the re- 
alization that shocks are likely to arise, at the 
latest when the fireball runs into an external 
medium, which would occur after the fireball is 
optically thin and would reconvert the kinetic en- 
ergy into nonthermal radiation. The complicated 
light curves can be understood in terms of in- 
ternal shocks in the outflow itself, before it 
runs into the external medium, caused by velocity 
variations in the outflow from the source. 

The next major breakthrough came in early 
1997 when the Italian-Dutch satellite Beppo-SAX 
succeeded in providing accurate X-ray measure- 
ments which, after a delay of 4-6 hours for pro- 
cessing, led to positions Q, allowing follow-ups 
at optical and other wavelengths, e.g. [127|. This 
paved the way for the measurement of redshift 
distances, the identification of candidate host 
galaxies, and the confirmation that they were in- 
deed at cosmological distances [^,^,|0[ . The de- 
tection of other GRB afterglows followed in rapid 
succession, sometimes extending to radio 24]25|] 
and over timescales of many months [12S], and in 
a number of cases resulted in the identification of 
candidate host galaxies, e.g. [1O8;D|,0, etc. The 
study of afterglows has provided strong confirma- 
tion for the generic fireball shock model of GRB. 
This model in fact led to a correct prediction j6l) , 
in advance of the observations, of the quantitative 
nature of afterglows at wavelengths longer than 
7-rays, which were in substantial agreement with 
the data ]l29| , p^ , pl^ , pm| , p^ . 

A major issue raised by the measurement of 
large redshifts, e.g. pO| , ^ , is that the mea- 
sured 7-ray fiuences imply a total energy of or- 
der 10^'^{ily /An) ergs, where Afl^ is the solid an- 
gle into which the gamma-rays are beamed. A 
beamed jet would clearly alleviate the energy re- 
quirements, but it is only recently that tentative 
evidence has been reported for evidence of a jet 
]5l| , p7| , p^ . Whether a jet is present or not, such 
energies are possible in the context of com- 
pact mergers involving neutron star-neutron star 
(NS-NS ) or black hole-neutron star (BH-NS ) bi- 
naries, or in hypernova/coUapsar models involv- 
ing a massive stellar progenitor |8iy95f| . In both 
cases, one is led to rely on MHD extraction of the 



spin energy of a disrupted torus and/or a central 
fast spinning BH, which can power a relativistic 
fireball resulting in the observed radiation. 

While it is at present unclear which, if any, 
of these progenitors is responsible for GRB, or 
whether perhaps different progenitors represent 
different subclasses of GRB, there is general 
agreement that they all would be expected to lead 
to the generic fireball shock scenario mentioned 
above. Much of the current effort is dedicated to 
understanding the progenitors more specifically, 
and trying to determine what effect, if any, they 
have on the observable burst and afterglow. 

2. Black Hole/Debris systems as generic 
GRB energy sources 

It has become increasingly apparent in the last 
few years that most plausible GRB progenitors 
suggested so far are expected to lead to a system 
with a central BH plus a temporary debris torus 
around it. Scenarios leading to this include, e.g. 
NS-NS or NS-BH mergers. Helium core - black 
hole [He/BH] or white dwarf - black hole [WD- 
BH] mergers, and the wide category labeled as hy- 
pernova or coUapsars including failed supernova 
lb [SNe lb], single or binary Wolf-Rayet [WR] 
collapse, etc. [pO 137.6^, |8l| , |95t , and accretion- 



induced collapse [L30,^. An important point 
is that the overall energetics from these various 
progenitors do not differ by more than about one 
order of magnitude |66[ . Another possibility is 
massive black holes (~ 10^ — 10^ Mq) in the ha- 
los of galaxies. Some related models involve a 
compact binary or a temporarily rotationally sta- 
bilized neutron star, perhaps with a superstrong 



field, e.g. [ [r25|p3|p3C|pI^ , which ultimately also 
should lead to a BH plus debris torus. 

Two large reservoirs of energy are available in 
these systems: the binding energy of the orbiting 
debris, and the spin energy of the black hole js^. 
The first can provide up to 42% of the rest mass 
energy of the disk, for a maximally rotating black 
hole, while the second can provide up to 29% of 
the rest mass of the black hole itself. The question 
is how to extract this energy. 

One energy extraction mechanisms is the vD —^ 
e+e~ process pO|, which can tap the thermal 
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energy of the torus produced by viscous dissi- 
pation. To be efRcient, the neutrinos must es- 
cape before being advected into the hole; on the 
other hand, the efficiency of conversion into pairs 
(which scales with the square of the neutrino 
density) is low if the neutrino production is too 
gradual. Typical estimates suggest a fireball of 



<10^^ erg |106,28|,^3[, except perhaps in the col- 
lapsar case where estimate 10^^'^ ergs for op- 
timum parameters. If the fireball is coUimated 
into a solid angle flj then of course the apparent 
"isotropized" energy would be larger by a fac- 
tor (47r/rjj) , but unless fij is <10~^ - IQ-* this 
would fail to satisfy the apparent isotropized en- 
ergy of 4 X lO^^^ ergs deduced for GRB 990123 

An alternative, and more efficient mechanism 
for tapping the energy of the torus may be 
through dissipation of magnetic fields generated 
by the differential rotation in the torus | |80| , |7l| , |62| , 
p3[ . Even before the BH forms, a NS-NS merging 
system might lead to winding up of the fields and 
dissipation in the last stages before the merger 
|55|Jl29|]. 

However, a larger energy source is available in 
the hole itself, especially if formed from a coalesc- 
ing compact binary, since then it is guaranteed to 
be rapidly spinning. Being more massive, it could 
contain more energy than the torus. The energy 
extractable in principle through MHD coupling to 
the rotation of the hole by the B-Z (Blandford & 
Znajek Q) effect could then be even larger than 

"ISlll. Col- 



that contained in the orbiting debris 
lectively, any such MHD outflows have been re- 
ferred to as Poynting jets. 

The various progenitors differ only slightly in 
the mass of the BH and that of the debris torus 
they produce, but they may differ more markedly 
in the amount of rotational energy contained in 
the BH. Strong magnetic fields, of order 10^^ 
G, are needed needed to carry away the rota- 
tional or gravitational energy in a time scale of 
tens of seconds [125,123|, which may be gener- 
ated on such timescales by a convective dynamo 
mechanism, the conditions for which are satis- 
fied in freshly collapsed neutron stars or neutron 
star tori If the magnetic fields do not 



thread the BH, then a Poynting outflow can at 
most carry the gravitational binding energy of 
the torus. For a maximally rotating and for a 
slow-rotating BH this is 

P - .AA .2 / 0.42(Mrf/MQ) ergs , (fast rot.); 
* - e^"oc I Q o6(Mrf/M0) ergs , (slow rot.), ' ^ > 

where e is the efficiency in converting gravita- 
tional into MHD jet energy. Th e torus or disk 
mass in a NS-NS merger is||l07| Md 10"^ - 



lO^-'Afo , and for a NS-BH, a He-BH, WD-BH 
merger or a binary WR collapse it may be esti- 
mated at [|l|,|§ Md IMq. In the He WD-BH 
merger and WR collapse the mass of the disk is 
uncertain due to lack of calculations on contin- 
ued accretion from the envelope, so IMq is just a 
rough estimate. The maximum torus-based MHD 
energy extraction is then 



8 X 10^^€{Md/MQ) ergs , 
1.2 X 10^^e{Md/MQ) ergs , 
0.8 X 10^^e{Md/0.1MQ) ergs 



for the NS-BH, He/WD-BH or collapsar case; the 
(slow rotating) failed SN lb case; and NS-NS case, 
respectively. 

If the magnetic fields in the torus thread the 
BH, the spin energy of the BH can in principle 
be extracted via the g (B-Z) mechanism (|2|). 
The extractable energy is 



Ebh - ef{a)AhhC^ 



(2) 



where e is the MHD efficiency factor, /(a) 
!-([! + Vl - a2]/2)i/2 < 0.29 is the rotational 
efficiency factor, and a ~ Jc/GM^ is the rota- 
tion parameter, which equals 1 for a maximally 
rotating black hole. The /(a) rotational factor 
is is small unless a is close to 1, where it rises 
sharply to its maximum value /(I) — 0.29, so 
the main requirement is a rapidly rotating black 
hole, a>,0.5. For a maximally rotating BH, the 
extractable energy is therefore 

E,naxM - 0.29eMbhC^ ~ 5xl0^^e{Mbh/MQ) ergs.(3) 

Rapid rotation is guaranteed in a NS-NS merger, 
since the radius (especially for a soft equation 
of state) is close to that of a black hole and 
the final orbital spin period is close to the re- 
quired maximal spin rotation period. The central 
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BH will have a mass 1 106, 107 1 of about 2.5Mq, 
so the NS-NS system can power a jet of up to 
Ens-ns <l.3xlO^'^e{Mbh/2.5MQ) ergs. A max- 
imal rotation rate may also be possible in a He- 
BH merger, depending on what fraction of the 
He core gets accreted along the rotation axis as 
opposed to along the equator |2^, and the same 
should apply to the binary fast-rotating WR sce- 
nario, which probably does not differ much in its 
final details from the He-BH merger. For a fast 
rotating BH of 2.5 — 3Mq threaded by the mag- 
netic field, the maximal energy carried out by the 
jet is then similar or somewhat larger than in 
the NS-NS case. The scenarios less likely to pro- 
duce a fast rotating BH are the NS-BH merger 
(where the rotation parameter could be limited 
to a < Mjis/Mbh, unless the BH is already fast- 
rotating) and the failed SNe lb (where the last 
material to fall in would have maximum angular 
momentum, but the material that was initially 
close to the hole has less angular momentum). 
The electromagnetic energy extraction from the 
BH in these could be limited by the /(a) fac- 
tor, but a lower limit would be given by the en- 
ergy available from the gravitational energy of the 
disk, in the second line of equation (0). 



log energy density (erg/c) 



Figure 1. Jet formation in a coUapsar model lead- 
ing to black hole from collapse of a fast rotating 
He core 15^1 



The total energetics differ thus between the 
various models at most by a factor 20 for a Poynt- 
ing (MHD) jet powered by the torus binding en- 
ergy, whereas for Poynting jets powered by the 
BH spin energy they differ at most by a factor of 
a few, depending on the rotation parameter. For 
instance, allowing for a total efficiency of 50%, a 
NS-NS merger whose jet is powered by the torus 
binding energy would require a beaming of the 7- 
rays by a factor {ATr/flj) ~ 100, or beaming by a 
factor 10 if the jet is powered by the B-Z mech- 
anism, to produce the equivalent of an isotropic 
energy of 4 x 10^'' ergs. The beaming require- 
ments of BH-NS and some of the collapsar scenar- 
ios are even less constraining, either when tapping 
the torus or the BH. Thus, even the most extreme 
energy requirements inferred observationally thus 
far can be plausibly satisfied by scenarios leading 
to a BH plus torus system. 

The major difference between the various mod- 
els is expected to be in the location where the 
burst occurs relative to the host galaxy (see 
They are also likely to differ substantially in the 
efficiency of producing a directly observable rela- 
tivistic outflow, as well as in the amount of col- 
limation of the jet they produce. The conditions 
for the efficient escape of a high-F jet are less pro- 
pitious if the "engine" is surrounded by an exten- 
sive envelope. In this case the jet has to "punch 
through" the envelope, and its ability to do so 
may be crucially dependent on the level of viscos- 
ity achieved in the debris torus (e.g. Q), higher 
viscosities leading to more powerful jets (see Fig- 
ure The simulations, so far, are nonrelativistic 
and one can only infer that high enough viscosi- 
ties can lead to jets capable of punching though 
a massive (several Mq) envelope. This is facili- 
tated, of course, if the envelope is fast-rotating, 
as in this case there is a centrifugally induced col- 
umn density minimum along the spin axis, which 
might be small enough to allow punch-through 
to occur. If they do, a very tightly coUimated 
beam may arise. "Cleaner" environments, such 
as NS-BH or NS-NS merger, or rotational sup- 
port loss / accretion induced collapse to BH would 
have much less material to be pushed out of the 
way by a jet, while their energy is, to order of 
magnitude, similar to that in massive stellar pro- 
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genitor cases. In these cases, on the other hand, 
there is no natural choke to coUimate a jet, which 
might therefore be somewhat wider than in mas- 
sive progenitor cases. 

3. The Fireball Shock Scenario 

Irrespective of the details of the progenitor, the 
resulting fireball is expected to be initially highly 
optically thick. From causality considerations the 
initial dimensions must be of order cit,ar <,10^ cm, 
where tyar is the variability timescale, and the 
luminosities must be much higher than a solar 
Eddington limit. Since most of the spectral en- 
ergy is observed above 0.5 MeV, the optical depth 
against 77 ^ is large, and an e^,7 fireball 
is expected. Due to the highly super-Eddington 
luminosity, this fireball must expand. Since in 
many bursts one observes a large fraction of the 
total energy at photon energies €y>lGeV, some- 
how the flow must be able to avoid degrading 
these photons (77 e"^ would lead, in a sta- 
tionary or slowly expanding flow, to photons just 
below 0.511 MeVQ ). In order to avoid this, 
it seems inescapable that the flow must be ex- 
panding with a very high Lorentz factor, since 
in this case the relative angle at which the pho- 
tons collide is less than and the threshold for 
the pair production is effectively diminished. The 
bulk Lorentz factor must be 



r^lO^(e7,ioGovet,Mcv)^^^ 



(4) 



in order for photons with energy e-y>,10 GeV to 
escape annihilation against target photons of en- 
ergy et ~ 1 MeV 34 1 . Thus, simply from 
observations and general physical considerations, 
a relativistically expanding fireball is expected. 
^From general considerations Q, one can see 
that an outfiow arising from an initial energy Eg 
imparted to a mass Mo << Eq/c^ within a ra- 
dius ri will lead to an expansion. Initially the 
bulk Lorentz F cx r, while comoving temperature 
drops cx r^^; however, F cannot increase beyond 

max 

^ ri Eo/MoC^, which is achieved at a ra- 
dius r/ri ~ 77, beyond which the flow continues 
to coast with F ^ 77 ~ constant [p8| . 



F - 



for r/ri<ri] 
for r/ri>^'q. 



However, the observed 7-ray spectrum observed 
is generally a broken power law, i.e., highly non- 
thermal. The optically thick e*7 fireball cannot, 
by itself, produce such a spectrum (it would tend 
rather to produce a modified blackbody, ||78| , ^ ). 
In addition, the expansion would lead to a con- 
version of internal energy into kinetic energy of 
expansion, so even after the fireball becomes op- 
tically thin, it would be highly inefficient, most 
of the energy being in the kinetic energy of the 
associated protons, rather than in photons. 

The most likely way to achieve a nonthermal 
spectrum in an energetically efficient manner is 
if the kinetic energy of the flow is re-converted 
into random energy via shocks, after the flow has 
become optically thin . This is a plausible sce- 
nario, in which two cases can be distinguished. In 
the first case (a) the expanding fireball runs into 
an external medium (the ISM, or a pre-ejected 
stellar wind [p8|j57|,^ , |109( |. The second possibil- 
ity (b) is that |9S|J77|| , even before external shocks 
occur, internal shocks develop in the relativistic 
wind itself, faster portions of the fiow catching 
up with the slower portions. This is a completely 
generic model, which is independent of the spe- 
cific nature of the progenitor, as long as it delivers 
the appropriate amount of energy (>,10^^ erg) in 
a small enough region (<,10'' cm). This model has 
been successful in explaining the major observa- 
tional properties of the gamma-ray emission, and 
is the main paradigm used for interpreting the 
GRB observations. 

External shocks will occur in an impulsive out- 
fiow of total energy Eo in an external medium of 
average particle density tIq at a radius and on a 
timescale 



tde 



^53 



rdec/icT') ~ 3 X lO^K 



-1/3 -2/3 
r?2 cm 

2^1/3^-1/3 



3/3 



(5) 



where the lab- frame energy of the swept-up exter- 
nal matter (F^WpC^ per proton) equals the initial 
energy Eg of the fireball, and 77 = F = 10^7^2 is 
the final bulk Lorentz factor of the ejecta. The 
typical observer-frame dynamic time of the shock 
(assuming the cooling time is shorter than this) 
is tdec ^ rdec/cT'^ ~ sccouds, for typical param- 
eters, and if, ~ tdec would be the burst duration 
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(the impulsive assumption requires that the ini- 
tial energy input occur in a time shorter than 
tdyn)- Variability on timescales shorter than tdec 
may occur on the cooling timescale or on the dy- 
namic timescale for inhomogeneities in the exter- 
nal medium, but generally this is not ideal for re- 
producing highly variable profiles[lll|. However, 
it can reproduce bursts with several peaks Q and 
may therefore be applicable to the class of long, 
smooth bursts. 

The same behavior T on r with conioving tem- 
perature cx , followed by saturation T^ax V 
at the same radius r/r; ~ rj occurs in a wind 
scenario ||79j , if one assumes that a lab-frame lu- 
minosity Lo and mass outflow Mq are injected at 
r ri and continuously maintained over a time 
tuj; here rj = Lo/MoC^- In such wind model, in- 
ternal shocks will occur at a radius and over a 
timescale lOSll 



^^tvar'n'^ ^ 3 X 10^'^t'(^arV2 Cm, 
ivar ^ ^dis / ) 



where shells of different energies Arj ~ rj initially 
separated by ctv (where ty < is the timescale 
of typical variations in the energy at r;) catch 
up with each other. In order for internal shocks 
to occur above the wind photosphere Vph ~ 
AfcrT/(47rTOpcr2) = 1.2 X lO^'^issry^^ cm, but 
also at radii greater than the saturation radius 
(so that most of the energy does not come out in 
the photospheric quasi-thermal radiation compo- 
nent) one needs to have 7.5 x 10^Lg(^i^a'/^ <,?73 x 

lO^Lgg^t^a/^. This type of models have the ad- 
vantage that they allow an arbitrarily compli- 
cated light curve, the shortest variation timescale 
i«or^lO~^ s being limited only by the dynamic 
timescale at r;, where the energy input may be 
expected to vary chaotically. Such internal shocks 
have been shown explicitly to reproduce (and be 
required by) some of the more complicated light 
curves 
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47,S7| (see however 



4. The Simple Standard Afterglow Model 

The dynamics of GRB and their afterglows can 
be understood in a fairly simple manner, inde- 
pendently of any uncertainties about the progen- 



itor systems, using a generalization of the method 
used to model supernova remnants. The simplest 
hypothesis is that the afterglow is due to a rel- 
ativistic expanding blast wave, which decelerates 
as time goes on |61 . The complex time structure 



of some bursts suggests that the central trigger 
may continue for up to 100 seconds, the 7-rays 
possibly being due to internal shocks. However, 
at much later times all memory of the initial time 
structure would be lost: essentially all that mat- 
ters is how much energy and momentum has been 
injected; the injection can be regarded as instan- 
taneous in the context of the much longer after- 
glow. As pointed in the original fireball shock pa- 
per [ p8[ , the external shock bolometric luminosity 
builds up and decays as 



L cx 



rise 



t-(i+9) decay 



(6) 



The first line is obtained equating, in the contact 
discontinuity frame, the kinetic flux L/47rr^ to 
the external ram pressure Pext^^ during the initial 
phase where T ~ constant, r cx i, while the second 
follows from energy conservation L (x E/t under 
adiabatic conditions (g takes into account radia- 
tive effects or bolometric corrections; the flux per 
unit frequency rises in the same way, and decays 
with (7 > 1 in equ. (||)). At the deceleration ra- 
dius (||) the fireball energy and the bulk Lorcntz 
factor decrease by a factor ~ 2 over a timescale 
tdec ~ i^dec/ (cT'^), and thereafter the bulk Lorentz 
factor decreases as a power law in radius. This is 

r cx r-ff cx t-f/(i+29) , r cx iVd+^s)^ (7) 

with g = (3,3/2) for the radiative (adiabatic) 
regime, in which pr^T ^ constant {pr^T^ ^ con- 
stant). At late times, a similarity solution [|| so- 
lution with 5 = 7/2 may be reached. 

The spectrum of radiation is likely to be due 
to synchrotron radiation, whose peak frequency 
in the observer frame is fm oc TB'j^, and both 
the comoving field B' and electron Lorentz fac- 
tor 7 are likely to be proportional to T . This 
implies that as T decreases, so will Vrn, and the ra- 
diation will move to longer wavelengths. This led 
1^,^ to early discussions, based on the forward 
blast wave, of the possibility of detecting at late 
times a radio or optical afterglow of the GRB. A 
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more detailed treatment of the fireball dynamics 
indicate that approximately equal amounts of en- 
ergy are radiated by the forward blast wave, mov- 
ing with ~ r into the surrounding medium, and 
by a reverse shock propagating with — 1 1 
back into the ejecta The electrons are there- 
fore shocked to much higher energies in the for- 
ward shock than in the reverse shock, producing a 
two-step synchrotron spectrum which during the 
deceleration time tdec peaks in the optical (re- 
verse) and in the "f/X (forward) p9| , |60[ . The pre- 
dicted fluences in the optical for typical bursts at 
cosmological distances were ~ 10~^'^ erg cm~^ 
s~^, or about a 9th magnitude prompt optical 
flash of duration comparable to the 7-rays, 
in agreement with a recent prompt optical de- 
tection in GRB 990123 §]. Detailed calculations 
and predictions of the time evolution of such a 
forward and reverse shock afterglow model ([|6l|) 
preceded the ob serva tions of the first afterglow 
GRB970228 (|P,|l27||), which was detected in 7- 
rays. X-rays and several optical bands, and was 
followed up for a number of months. 

The simplest spherical afterglow model gener- 
ally concentrates only on the properties of the 
forward blast wave radiation, for which the flux 
at a given frequency and the synchrotron peak 
frequency decay at a rate ]6l| , |67|| 

CX f[3-2ff(l-2/3)]/(l+2<,) ^ ^^^^-4W(l+2<,)^(8) 

where g is the exponent of F (equ. 0) and /3 is the 
photon spectral energy slope. The decay rate of 
the forward shock Fi, in equ.(^) is typically slower 
than that of the reverse shock , and the rea- 
son why the " simplest" model was stripped down 
to its forward shock component only is that, for 
the first two years 1997-1998, afterglows were fol- 
lowed in more detail only after the several hours 
needed by Beppo-SAX to acquire accurate posi- 
tions, by which time both reverse external shock 
and internal shock components are expected to 
have become unobservable. This simple standard 
model has been remarkably successful at explain- 
ing the gross features and light curves of GRB 
970228, GRB 970508 (after 2 days; for early rise, 
see S||) e.g. [ p6|pTlJl3^JT0l|] (see Figure |). 

This simplest afterglow model produces at any 
given time a three-segment power law spectrum 
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Figu re 2. GRB 970228 hght-curves compared 
|13£] to the blast wave model predictions of ||6l|] 



with two breaks. At low frequencies there is 
a steeply rising synchrotron self-absorbed spec- 
trum up to a self-absorption break Va, followed 
by a -1-1/3 energy index spectrum up to the syn- 
chrotron break corresponding to the minimum 
energy 7^ of the power-law accelerated electrons, 
and then a —{p— l)/2 energy spectrum above 
this break, for electrons in the adiabatic regime 
(where 7"^ is the electron energy distribution 
above 7m). A fourth segment and a third break 
is expected at energies where the electron cool- 
ing time becomes short compared to the expan- 
sion time, with a spectral slope — p/2 above that. 
With this third "cooling" break i^b, first calcu- 
lated in 1^ and more explicitly detailed in [112|, 
one has what has come to be called the simple 
"standard" model of GRB afterglows. One of the 
predictions of this model is that the relation 
between the temporal decay index a, for g = 3/2 
in F cx r~3, is related to the photon spectral en- 
ergy index (3 through through 



cx Tj/^ ,with a = (3/2)/3 



(9) 



This relationship appears to be valid in many (al- 
though not all) cases, especially after the first few 
days, and is compatible with an electron spectral 
index p ~ 2.2 — 2.5 which is typical of shock accel- 
eration, e.g. |132,112,135|, etc. As the remnant 
expands the photon spectrum moves to lower fre- 
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quencies, and the flux in a given band decays as 
a power law in time, whose index can change as 
breaks move through it. For the simple standard 
model, snapshot overall spectra have been de- 
duced by extrapolating spectra at different wave- 
bands and times using assumed simple time de- 
pendences 1 133, 135 1 . These can be used to derive 
rough fits for the different physical parameters of 
the burst and environment, e.g. the total energy 
E, the magnetic and electron-proton coupling pa- 
rameters cb and and the external density rio 
(see Figure p|) . 
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Figure 3. Snapshot spectrum of GRB 970508 at 
t = 12 days and standard afterglow model fit |135] 



Since the simple afterglow model has generally 
proved quite robust, it is worth reviewing the 
assumptions made in it. The following apply 
both to the "simple standard" model using for- 
ward shocks only, and to the original |6l| version 
including both forward and reverse shocks: 

a) A single value of Eo and To — rj is used, 

b) the external medium riext is homogeneous, 

c) the accelerated electron spectral index p, the 
magnetic field and electron to proton equipartion 
ratios Sb and Eg do not change in time, 

d) the expansion is relativistic and the dynamics 
are given by F cx r~^/^ (adiabatic), 

d) the outflow is spherical (or angle independent 
inside some jet solid angle fij), 

e) the observed radiation is characterized by the 
scaling relations along the line of sight. 



These assumptions, even if correct over some 
range, clearly would break down after some 
time. Estimates for the time needed to reach 
the non-relativistic expansion regime are typi- 
cally <, month(s) ( |l2g| ]), or less if there is an 
initial radiative regime F cx r~^. However, even 
when electron radiative times are shorter than the 
expansion time, it is unclear whether a regime 
F cx should occur, since it would require 
strong electron-proton coupling As far as 

sphericity, the standard model can be straight- 
forwardly generalized to the case where the en- 
ergy is assumed to be channeled initially into a 
sohd angle fl^ < An In this case |T04|jlO^ ] 

a change occurs after F drops below il. 



-1/2 



af- 



ter which the side of the jet becomes observable, 
and soon thereafter one expects a faster decay 
of F if the jet starts to expands sideways, lead- 
ing to a decrease in the brightness. A calculation 
based on the sideways expansion, using the usual 
scaling laws for a single central line of sight |105] 
leads then to a steepening of the light curve. Un- 
til recently, no evidence for a steepening could be 
found in afterglows over several months. E.g., in 
GRB 971214 a snapshot standard model fit 
and the lack of a break in the late light curve 
could be, in principle, interpreted as evidence for 
lack of a jet, leading to an (isotropic) energy esti- 
mate of 1053-5 ergs. While such large energy out- 
puts are possible in either NS-NS, NS-BH merg- 
ers or in hypernova/coUapsar models |81 9^ 
using MHD extraction of the spin energy of a dis- 
rupted torus and/or a central fast spinning BH, 
it is worth stressing that what these snapshot fits 
constrain is only the energy per solid angle ]6^ . 
Also, the expectation of a break after some weeks 
or months (e.g., due to F dropping either below a 



few, or below Q 



-l/2^ 



is based upon the simple im- 



pulsive (angle-independent delta or top-hat func- 
tion) energy input approximation. The latter is 
useful, but departures from it would be natural, 
and certainly not surprising. In fact, as discussed 
below, tentative evidence for beaming in one obe- 
jct has recently been reported pl| , p7|JT2| , but it is 
difficulty to ascertain, and could be masked by a 
number of commonly expected effects. 
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Figure 4. Ring-like equal-arrival time T surfaces 
of an afterglow, based on pa] 



5. "Post-standard" Afterglow Models 

In a realistic situation, one could expect any of 
several fairly natural departures from the simple 
standard model to occur. The first one is that 
the emitting region seen by the observer resem- 
bles a ring [p^p||lTo| (see Figure |). This ef- 
fect may, in fact, be important in giving rise to 
the radio scintillation pattern seen in several af- 
terglows, since this requires the emitting source 
to be of small dimensions, which is aided if the 
emissio n is r ing-like, e.g. in the example of GRB 
970508 @ (Figure I). 




Figure 5. Radio afterglow light-curves of 
GRB970508 at 4.86GHz and 1.43GHz, compared 
with the predictions of [134|. 



One expects afterglows to show a diversity in 
their decay rates, not only due to different /3 but 



also from the possibility of a non-standard rela- 
tion between the temporal decay index a and the 
spectral energy index /?, different from equ. (||). 
The most obvious departure from the simplest 
standard model occurs if the external medium is 
inhomogeneous: for instance, for n^^t oc r"'', the 
energy conservation condition is Y'^r'~'^ ^ con- 
stant, which changes significantly the temporal 
decay rates |^]. Such a power law dependence 
is expected if the external medium is a wind, say 
from an evolved progenitor star as implied in the 
hypernova scenario (such winds are generally used 
to fit supernova remnant models). Another obvi- 
ous non-standard effect, which it is reasonable to 
expect, is departures from a simple impulsive in- 
jection approximation (i.e. from a delta or top 
hat function with a single value for Eo and Fq). 
An example is if the mass and energy injected 
during the burst duration (say tens of sec- 
onds) obeys M(> F) cx F"", E(> F) cx F^-", i.e. 
more energy emitted with lower Forentz factors at 
later times (but still shorter than the gamma-ray 
pulse duration). This would drastically change 
the temporal decay rate and extend the after- 
glow lifetime in the relativistic regime, providing 
a late "energy refreshment" to the blast wave on 
time scales comparable to the afterglow time scale 
I IOC]. These two cases lead to a decay rate 



F cx r ^ 



j.-(3-<i)/2 
y.-3/(2+s) 



; next (xr 

; E{> F) cx r^' 



(10) 



Expressions for the temporal decay index 
a(/3, s, d) in F^, cx are given by [6^JlO(]| ], which 
now depend also on s and/or d (and not just on /3 
as in the simple standard relation of equ.(||). The 
result is that the decay can be flatter (or steeper, 
depending on s and d) than the simple standard 
a — (3/2)/3. A third non-standard effect, which is 
entirely natural, occurs when the energy and/or 
the bulk Lorentz factor injected are some function 
of the angle. A simple case is Eg cx 6^^ , Fq cx 9^'^ 
within a range of angles; this leads to the outflow 
at different angles shocking at different radii and 
its radiation arriving at the observed at differ- 
ent delayed times, and it has a marked effect on 
the time dependence of the afterglow ]^, with 
a — a{/3,j,k) flatter or steeper than the stan- 
dard value, depending on j, k. Thus in general, a 
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Figure 6. Optical light-curve of GRB 970508, fit- 
ted with a non-uniform injection model (a simi- 
lar fit can be obtained with an off-axis jet plus a 
weaker isotropic component) ]8q] 



temporal decay index which is a function of more 
than one parameter 

F^oif^v^ ,with a = a(/3, d, s, j, fc, • • •) , (11) 

is not surprising; what is more remarkable is that, 
in many cases, the simple standard relation (||) is 
sufficient to describe the gross overall behavior at 
late times. 

Strong evidence for departures from the simple 
standard model is provided by, e.g., sharp rises 
or humps in the light curves followed by a re- 
newed decay, as in GRB 970508 (J8|,Q). De- 
tailed time-dependent model fits | |88| to the X- 
ray, optical and radio light curves of GRB 970228 
and GRB 970508 show that, in order to explain 
the humps, a non-uniform injection (Figure ^) or 
an anisotropic outflow is required. 

These fits indicate that the shock physics may 
be a function of the shock strength (e.g. the elec- 
tron index p, injection fraction C, and/or e;,, eg 
change in time), and also indicate that dust ab- 
sorption is needed to simultaneously fit the X-ray 
and optical fluxes. The effects of beaming (out- 
flow within a limited range of solid angles) can 
be significant |Q , but are coupled with other ef- 
fects, and a careful analysis is needed to disen- 
tangle them. 

One consequence of "post-standard" decay 



laws (e.g. from density inhomogeneities, non- 
uniform injection or anisotropic outflow) is that 
the transition to a steeper jet regime F < ~ 
few can occur as late as six months to a year after 
the outburst, depending on details of the energy 
input. This transition is made more difficult to 
detect by the fact that, as numerical integration 
over angles of the ring-like emission show, the 
transition is very gradual and the effects of side- 
ways expansion effects are not so drastic as in- 
ferred [105 1 from the scaling laws along the central 
line of sight. This is because even though the flux 
from the head-on part of the remnant decreases 
faster, this is more than compensated by the in- 
creased emission measure from sweeping up ex- 
ternal matter over a larger angle, and by the fact 
that the extra radiation, arising at larger angles, 
arrives later and re-fills the steeper light curve. 
The inference (e.g. [p7p05|) that GRB 970508 
and a few other bursts were isotropic due to the 
lack of an observable break is predicated entirely 
on the validity of the simplest standard fireball 
assumption. Since these assumptions are dras- 
tic simplifications, and physically plausible gen- 
eralizations lead to different conclusions, one can 
interpret the results of [^,105| as arguments in- 
dicating that post- standard features are, in fact, 
necessary in some objects. 

6. Prompt multi-wavelength flashes, re- 
verse shocks and jets 

Prompt optical. X-ray and GeV flashes from 
reverse and forward shocks, as well as from in- 
ternal shocks, have been calculated in theoreti- 
cal fireball shock models for a number of years 

f, |60| , p2| , |6l|J113| ], as have bee n jets (e.g . |55 5^, 
, and in more detail |104, 8^ , |8^ , |l05| ). Thus, 



while in recent years they were not explicitly part 
of the "simple standard" model, they are not 
strictly "post" -standard either, since they gener- 
ally use the "standard" assumptions, and they 
have a long history. However, observational evi- 
dence for these effects were largely lacking, until 
the detection of a prompt (within 22 s) optical 
flash from GRB 990123 with ROTSE by §, to- 
gether with X-ray, optical and radio follow-ups 
citekul99,gal99,fru99,and99,cas99,hjo99. GRB 
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990123 is so far unique not only for its prompt 
optical detection, but also by the fact that if it 
were emitting isotropically, based on its redshift 
z = 1.6 its energy would be the largest of 

any GRB so far, 4 x 10^** ergs. It is, however, 
also the first (tentative) case in which there is 
evidence for jet-like emission ^1 27,1^]. An ad- 
ditional, uncommon feature is that a radio after- 
glow appeared after only one day, only to disap- 
pear the next |3^j5l[ |. 

The prompt optical hght curve of GRB 990123 
decays initially as oc t~^-^ to cx t"^-^ Q, much 
steeper than the typical cx of previous opti- 
cal afterglows detected after several hours. How- 
ever, after about 10 minutes its decay rate moder- 
ates, and appears to join smoothly onto a slower 
decay rate cx t~^-^ measured with large telescopes 
|30[|5|,|2|,|l|l after hours and days. The prompt 
optical flash peaked at 9-th magnitude after 55 s 
|2|, and in fact a 9-th magnitude prompt flash 
with a steeper decay rate had been predicted 
more than two years ago ||6l[| , from the syn- 
chrotron radiation of the reverse shock in GRB af- 



terglows at cosmological redshifts (see also |113]). 
An optical flash contemporaneous with the 7-ray 
burst, coming from the reverse shock and with flu- 
ence corresponding to that magnitude, had also 
been predicted earlier |5|,^^. An origin of the 
optical prompt flash in internal shocks |]6l| , |67| 
cannot be ruled out yet, but is less likely since 
the optical light curve and the 7-rays appear not 
to correlate well |114,Q (but the early optical 
light curve has only three points). The subse- 
quent slower decay agrees with the predictions 
for the forward component of the external shock 
6lUll^j6^. 

The evidence for a jet is possibly the most ex- 
citing, although must still be considered tenta- 
tive. It is based on an apparent steepening of 
the light curve after about three days | |5l| , |27| , p^ . 
This is harder to establish than the decay of the 
two previous earlier portions of the light curve, 
since by this time the flux has decreased to a 
level where the detector noise and the light of 
the host galaxy become important. However, af- 
ter correcting for this, the r-band data appears 
to steepen significantly. (In the K-band, where 
the noise level is higher, a steepening is not obvi- 



ous, but the issue should be settled with further 
Space Telescope observations). If real, this steep- 
ening is probably due to the transition between 
early relativistic expansion, when the light-cone 
is narrower than the jet opening, and the late ex- 
pansion, when the light-cone has become wider 
than the jet, leading to a drop in the effective 
flux ||T0^,|5l|,|67|,p^ . A rough estimate leads to 
a jet opening angle of 3-5 degrees, which would 
reduce the total energy requirements to about 
4 X 10^^ ergs. This is about two order of mag- 
nitude less than the binding energy of a few solar 
rest masses, which, even allowing for substantial 
inefficiencies, is compatible with currently favored 
scenarios (e.g. ^,Q) based on a stellar collapse 
or a compact binary merger. 

7. Location and Environmental Effects 

The location of the afterglow relative to the 
host galaxy center can provide clues both for the 
nature of the progenitor and for the external den- 
sity encountered by the fireball. A hypernova 
model would be expected to occur inside a galaxy 
in a high density environment Ug > 10^ — 10^ 
cm^"^. Most of the detected and well identified 
afterglows are inside the projected image of the 
host galaxy 0, and some also show evidence for 
a dense medium at least in front of the after- 
glow (1^^). For a number of bursts there are 
constraints from the lack of a detectable, even 



faint, host galaxy |116], but at least for Beppo- 
SAX bursts (which is sensitive only to long bursts 
i6>,20 s) the success rate in finding candidate 
hosts is high. 

In NS-NS mergers one would expect a BH plus 
debris torus system and roughly the same total 
energy as in a hypernova model, but the mean 
distance traveled from birth is of order several 
Kpc [|l0| , leading to a burst presumably in a less 
dense environment. The fits of 1 135] to the obser- 
vational data on GRB 970508 and GRB 971214 
in fact suggest external densities in the range of 
Uo ~ 0.04-0.4 cm~^, which would be more typical 
of a tenuous interstellar medium. These could be 
within the volume of the galaxy, but for NS-NS 
on average one would expect as many GRB inside 
as outside. This is based on an estimate of the 
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mean NS-NS merger time of 10® years; other esti- 



mated merger times (e.g. 10'' years, |l26| ) would 
give a burst much closer to the birth site. BH- 
NS mergers would also occur in timescales <,10^ 
years, and would be expected to give bursts well 
inside the host galaxy ([|lO|; see however |^). In 
at least one "snapshot" standard afterglow spec- 
tral fit for GRB 980329 [|lO| the deduced external 



density is Uo ~ 10 cm~'^. In some of the other de- 
tected afterglows there is other evidence for a rel- 
atively dense gaseous environments, as suggested, 
e.g. by evidence for dust |l02| in GRB970508, 
the absence of an optical afterglow and presence 
of strong soft X-ray absorption |^,|9| in GRB 
970828, the lack an an optical afterglow in the 



(radio-detected) afterglow (|22|) of GRB980329, 
and spectral fits to the low energy portion of the 
X-ray afterglow of several bursts The lat- 

ter observations may be suggestive of hypernova 
models , involving the collapse of a massive 

star or its merger with a compact companion. 

One important caveat is that all afterglows 
found so far are based on Beppo-SAX positions, 
which is sensitive only to long bursts t6>,20 s |39| . 
This is significant, since it appears likely that NS- 
NS mergers lead Q to short bursts with tb<lO 
s. To make sure that a population of short GRB 
afterglows is not being missed will probably need 
to await results from HETE iBSl and from the 



planned Swift |12C] mission, which is designed to 
accurately locate 300 GRB/yr. 

An interesting case is the apparent coincidence 
of GRB 980425 with the unusual SN Ib/Ic 1998bw 
pgj , which may represent a new class of SN [po|p| . 
If true, this could imply that some or perhaps all 
GRB could be associated with SN Ib/Ic ||l3^, dif- 
ferring only in their viewing angles relative to a 
very narrow jet. Alternatively, the GRB could be 
(e.g. 1 138 1) a new subclass of GRB with lower 
energy ~ {flj / Att) erg, only rarely ob- 
servable, while the great majority of the observed 
GRB would have the energies ~ 10^''(Slj/47r) 
ergs as inferred from high redshift observations. 
The difhculties are that it would require extreme 
collimations by factors 10~^ — lO""*, and the sta- 
tistical association is so far not significant Q]. 

The environment in which a GRB occurs 
should also influence the nature of the afterglows 



in other ways. The blast wave and reverse shock 
that give rise to the X-rays, optical, etc occur over 
timescales proportional to tdec « n^^/^ (equ.Q) 
which is longer in lower density environments, so 
for the same energy the flux is lower, roughly 

1/2 

Fi^ cx EoTiJ^^ , contributing also to make after- 
glows in the intergalactic medium harder to de- 
tect. However, in addition to affecting broad- 
band fluxes, one may also expect specific spectral 
signatures from the external medium imprinted in 
the X-ray and optical continuum, such as atomic 
edges and lines [p pl|j6^ . These may be used 
both to diagnose the chemical abundances and 
the ionization state (or local separation from the 
burst), as well as serving as potential alternative 
redshift indicators. (In addition, the outflowing 
ejecta itself may also contribute blue-shifted edge 
and line features, especially if metal-rich blobs or 
filaments are entrained in the flow from the dis- 
rupted progenitor debris |]63|| , which could serve 
as diagnostic for the progenitor composition and 
outflow Lorentz factor). To distinguish between 



progenitors 
see also [ml 




, an interesting prediction ( |64| ; 
is that the presence of a measur- 
able Fe K-a X-ray emission line could be a diag- 
nostic of a hypernova, since in this case one may 
expect a massive envelope at a radius comparable 
to a light-day where tt<1^ capable of reprocess- 
ing the X-ray continuum by recombination and 
fluorescence. Two groups p4,139| have in fact re- 



cently reported the possible detection of Fe emis- 
sion fines in GRB 970508 and GRB 970828. 

8. Conclusions 

The fireball shock model of gamma-ray bursts 
has proved quite robust in providing a consis- 
tent overall interpretation of the major features 
of these objects at various frequencies and over 
timescales ranging from the short initial burst 
to afterglows extending over many months. The 
standard internal shock scenario is able to re- 
produce the properties of the 7-ray light curves, 
while external shocks involving a forward blast 
wave and a reverse shock are successful in re- 
producing the afterglows observed in X-rays, op- 
tical and radio. The "simple standard model" 
of afterglows, involving four spectral slopes and 



13 



three breaks is quite useful in understanding the 
'snapshot' multiwavelength spectra of most after- 
glows. However, the effects associated with a jet- 
like outflow and the possible differential beaming 
at various energies requires further investigations, 
both theoretical and observational. Caution is re- 
quired in interpreting the observations on the ba- 
sis of the simple standard model. For instance, 
more detailed numerical models, as opposed to 
the more common analytical scaling law mod- 
els, show that the contributions of radiation from 
different angles and the gradual transition be- 
tween different dynamical and radiative regimes 
lead to a considerable rounding-off of the spec- 
tral shoulders and light-curve slope changes, so 
that breaks cannot be easily located unless the 
spectral sampling is dense and continuous, both 
in frequency and in time. Some of the observed 
light curves with humps, e.g. in GRB 970508, re- 
quire 'post-standard' model features (i.e. beyond 
those assumed in the standard model), such as ei- 
ther non-uniform injection episodes or anisotropic 
outflows. Time-dependent multiwavelength fits 
jssl of some bursts also indicate that the parame- 
ters characterizing the shock physics change with 
time. Even without humps or slope changes, a 
non-standard relation between the spectral and 
temporal decay slope is observed in several ob- 
jects, e.g. in GRB 990123 |l[. These are, in our 
view 1^ , a strong indication for " post-standard" 
effects in such bursts. 

Much progress has been made in understand- 
ing how gamma-rays can arise in fireballs pro- 
duced by brief events depositing a large amount 
of energy in a small volume, and in deriving the 
generic properties of the long wavelength after- 
glows that follow from this. There still remain 
a number of mysteries, especially concerning the 
identity of their progenitors, the nature of the 
triggering mechanism, the transport of the en- 
ergy, the time scales involved, and the nature 
and effects of beaming. However, even if we do 
not yet understand the details of the gamma-ray 
burst central engine, it is clear that these phe- 
nomena are among the most powerful transients 
in the Universe, and they could serve as power- 
ful beacons for probing the high redshift {z > 5) 
universe. The modeling of the burst mechanism 



itself, as well as the resulting outflows and radia- 
tion, will continue to be a formidable challenge to 
theorists and to computational techniques. How- 
ever, the theoretical understanding appears to be 
converging, and with dedicated new and planned 
observational missions under way, the prospects 
for significant progress are realistic. 

I am grateful to Martin Rees for stimulating 
collaborations, as well as to Alin Panaitescu, Hara 
Papathanassiou and Ralph Wijers. 
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